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Abstract 

We analyze the excitation of Ar substrate in contact with Na clusters using a previously developed hierarchical 
model for the description of the system cnstituted of a highly reactive metal cluster in contact with a rather inert 
substrate. Particular attention is paid to the dipole excitation of the Ar atoms and the energy stored therein. The 
Na clusters are considered at different charge states, anions, cations, and neutral clusters for the case of deposition 
and a highly ionized cluster embedded in a matrix. It is found that the dipole polarization of the Ar atoms stores 
the largest fraction of energy in the case of charged clusters. Some, although smaller, polarization is also observed 
for polar clusters, as Nas. The effect is predominantly induced by the electrostatic interaction. 
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1. Introduction 

The study of clusters in contact with an envi- 
ronment has motivated many investigations over 
the years. Two major situations can be considered, 
namely a cluster deposited on a surface or a clus- 
ter embedded inside a matrix. Both situations bear 
some similarities and bring complementing infor- 
mation. Clusters on surfaces provide interesting 
perspectives for basic research and for applications 
to nano-structured materials |H2j . e.g., the synthe- 
sis of deposited clusters, either controlled growth of 
elementary units on a surface by molecular beam 
epitaxy [3] or direct deposition of size-selected clus- 
ters on a substrate [3]. One should also mention 
the non-destructive deposition technique of metal 
clusters on metal surfaces using a thin rare gas film 
above the metal surface [5] . At the side of embedded 
species, one can refer to the many studies exploring 
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the optical response in various dynamical scenarios 
|6|7|8j complementing similar investigations in the 
case of free clusters. Note furthermore that embed- 
ded clusters (or molecules) may also be viewed as 
model systems for a detailed analysis of radiation 
effects in matter |9ll0j . 

^From the theory side, the description of clus- 
ters in contact with an environment implies an 
extra complexity because one needs to account 
for the degrees of freedom of the substrate. The- 
oretical descriptions thus predominantly employ 
classical molecular dynamics with effective atom- 
atom forces, see [IT]. This was, for example, done 
for the deposition dynamics of Cu clusters on 
metal [T5] or Ar [T3] surfaces, and of Al or Au 
clusters on Si02 [Hj. Such simplified approaches 
overlook possible effects from electronic degrees of 
freedom which can become crucial in metal clus- 
ters, particularly if a finite net charge is involved. 
It is thus often necessary to account explicitely 
for electronic degrees of freedom. Fully detailed 
calculations have been performed which treat all 
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constituents with their electronic dynamics, e.g. 
for the structure of small Na clusters on NaCl [15J 
or the deposit dynamics of Pd clusters on a MgO 
substrate [16]. But the numerical effort of such 
detailed computations quickly grows huge. Further- 
more, these subtle models are hardly extendable 
to truly dynamical situations, to larger clusters or 
substrates, and to systematic explorations for broad 
variations of conditions. This leaves space for a 
manifold of approximations allowing an affordable 
compromise between reliability and expense. Such 
methods are often called quantum-mechanical- 
molecular- mechanical (QM/MM) models and have 
been applied for instance to chromophores in bio- 
molecules |17|18j . surface physics |19|20j . materials 
physics |21l22l23l24j . embedded molecules [25J and 
ion channels of cell membranes [25]. We have de- 
veloped such a QM/MM modeling, primarily in the 
case of Na in contact with Ar 27 28 29J and success- 
fully applied this method to deposition dynamics 
on finite Ar clusters [515] or on Ar surfaces [31], as 
well as to irradiation scenarios in the case of em- 
bedded clusters both in the linear (optical response 
|32j ) and non linear (hindered explosion [33134] ) 
domains. More recently, we have extended the mod- 
eling to a MgO substrate [35 . The originality of 
this approach lies in the fact that the environment 
polarizability is treated dynamically, a key aspect 
especially when charged species are considered [36 . 

In this paper, we want to further investigate the 
importance of this dynamical treatment of the envi- 
ronment's polarizability by exploring how the sub- 
strate's dipoles respond to a perturbation. We shall 
consider two typical scenarios : Cluster deposition 
on a surface and irradiation of an embedded cluster. 
In both cases, we shall explore in particular the im- 
pact of charge, either because the deposited species 
is charged, or because the embedded cluster acts as 
a chromophorc in a laser field and thus quickly ac- 
quires charge after irradiation. This exploration is 
focused on theoretical aspects but is qualitatively 
closely related to recent experiments in which the de- 
position of charged Ag clusters on a Ar matrix, itself 
deposited on a Au surface has been studied [37138] . 
It was shown in this experiment that the substrate, 
in spite of its a priori inert character, acquires a sub- 
stantial inner excitation which plays a key role in 
the whole process. Exploring the dipole degrees of 
freedom in our Na-Ar combination then represents 
the simplest model case in relation to these experi- 
ments. We shall see that indeed, as suggested by our 
earlier investigations [36], these internal degrees of 



freedom are readily excited in most of the studied 
cases. 

The paper is organized as follows. After a brief 
reminder of the content of the model used here, we 
quickly focus on relevant test cases. We first show 
the overall importance of charge effects and then 
study the spatial extension of the inner excitation of 
the environment. We take examples from deposition 
processes and irradiation of embedded clusters. 

2. Model 

We first give a very brief summary of the hierar- 
chical description of the combined Na-Ar system. 
We treat the metal atoms in full microscopic detail 
at the level of Time Dependent Local Density Ap- 
proximation (TDLDA) for the valence electrons of 
the Na cluster. An average self-interaction correc- 
tion is applied to LDA in order to put the ioniza- 
tion potential of the Na system at the correct place 
[35] . Electronic degrees of freedom are coupled to 
Molecular Dynamics (MD) for the ions. Details on 
this very successful TDLDA-MD approach for free 
clusters can be found in [40141] , The environment 
(substrate or matrix) consists out of Ar atoms which 
are described by classical degrees of freedom, both 
in terms of position and dipole moment. The latter 
serves to take into account the dynamical polariz- 
ability of the substrate atoms. They have been ad- 
justed carefully to recover the static and dynamical 
polarizabilities of Ar atoms. Electronic emission is 
not possible at the side of Ar atoms which limits the 
violence of the processes studied. But it is easy to 
check, in terms of the amplitude of Ar dipoles, how 
much energy is actually absorbed by an Ar atom. 
More precisely, this internal excitation energy of the 
Ar atoms is related, in our model, to the Ar dipole 
amplitude d by 

^di P = ^e 2 ^d 2 , (1) 

2 C*Ar 

where o^r is the static polarizability of bulk Ar, and 
QAr the effective charge of the Ar cores [25] . 

As long as this energy (and correspondigly the 
dipole amplitude) remains safely below the Ar atom 
ionization potential, the model is perfectly applica- 
ble. We checked that this is always the case in all sit- 
uations encountered here. For then, the dipole am- 
plitude will represent the only possible inner excita- 
tion of Ar atoms. The Ar atoms are coupled to the 
Na by a long range polarization potential and some 
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short range repulsion to account for the Pauli block- 
ing of cluster electrons in the vicinity of the Ar cores. 
The model has been calibrated to measured prop- 
erties of typical Na-Ar systems. We refer the reader 
to [27128142] for a detailed description of the model 
and the detailed fitting of the various parameters. 

In the case of deposition dynamics, the Ar(OOl) 
surface is modelled through an Ai'384 system, con- 
sisting in six layers of 8x8 Ar atoms with atoms in 
the two lowest layers frozen at bulk crystal positions. 
The layers are furthermore periodically repeated in 
both lateral directions, thus simulating bulk mate- 
rial in these two dimensions. We have checked that 
the finiteness of the sampling does not alter the re- 
sults, at least qualitatively 30 36J . The dynamics is 
then initialized by placing the projectile at a finite 
distance (15 ao) from the surface and boosting it 
with a given initial kinetic energy Eq, towards the 
substrate and along the direction normal to it (de- 
noted by z in the following). In the case of embed- 
ded species, the system is constructed by consider- 
ing first a finite piece of bulk Ar, drilling a small hole 
in the center to insert the Nas cluster and then op- 
timizing the whole ensemble (cluster + embedding 
material Ar434) at the side of cluster electrons and 
ions and in terms of Ar positions and dipoles. The 
system is then irradiated by a laser and its response 
followed in time. We analyze the subsequent dynam- 
ics in terms of detailed ionic and atomic coordinates 
and dipoles. 



3. Impact of charge 

Let us first consider a typical example, namely the 
deposition of a Na^~ cluster (consisting in a pentagon 
with a top ion on its symmetry axis) on Ar(001). 
We take a charged cluster in the spirit of [37138] 
and a moderate initial kinetic energy in order to 
observe sticking of the cluster on the surface with 
no destruction of the substrate |31|43j . Fig.[l]shows 
the time evolution of positions and energies. In the 
upper panel are plotted the z-coordinates of ions 
(Nag ) and surface (Ar) atoms. The z direction is 
the direction along the initial cluster velocity, per- 
pendicular to the surface. The Na^ cluster is sig- 
nificantly perturbed as one can see from the large 
oscillations of the initially top ion going through 
the pentagon plane. Nonetheless, the overall process 
converges steadily to a robust sticking of the cluster 
on the surface. The surface itself is perturbed with 
some ionic rearrangement but with preservation of 
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Fig. 1. Time evolution of coordinates and energies during dy- 
namical deposition of Nag" on Ar(001) with initial kinetic en- 
ergy Eg = 0.136 eV/ion. Upper panel: z coordinates. Lower 
panel: total kinetic energy of Ar cores, Na ions and total 
dipole excitation energy. 

layer structure. The lower panel of Fig. [T] displays 
the time evolution of energies, i.e. kinetic energies of 
cluster ions and Ar atoms, and the energy stored in 
Ar dipoles, see Eq. (JlJ . The time evolution is rather 
simple with an almost instantaneous transfer of clus- 
ter kinetic energy to substrate degrees of freedom 
30 31]. Note that some part of initial energy is flow- 
ing into potential energy which is not shown here. 
Still, the interesting feature is that, while Ar atoms 
acquire a significant kinetic energy of about half the 
maximum one of the cluster, they store about twice 
as much energy in their dipoles. This is a key as- 
pect. It means that the substrate is not only heated 
up by cluster impact but also internally excited at 
the side of each constituent atom. This shows that 
a proper treatment of the dynamical surface polar- 
izability can be crucial, particularly if a charge is 
involved. The energy sharing is established almost 
instantaneously at impact time and energies then 
remain rather constant in time. 

The strong effect at the side of dipoles observed 
in Fig. [I] has to be explored further in order to try 
to identify where it comes from. Using a Nag" clus- 
ter as a projectile implies two possible effects, from 
charge and/or mass. We have considered here the 
case of a cationic cluster but an anionic one could 
be envisioned as well. We thus try to disentangle 
charge and mass effects by considering various pos- 
sible combinations namely Nag , Na^ , Nag , Na, Na + 
and Na - . We hence focus on dipole energies of the 
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six test cases, which are plotted as a function of time 
in Fig. [5J As compared to the Na^ case, the dipole 
energies are suppressed by several orders of magni- 
tude (2 for Nag and even 6 for Na) for neutral projec- 
tiles, while with charged projectiles, they are of the 
same order, with even a higher dipole energy for a 
positive charge by a factor about 1.5 with respect to 
the Ar core kinetic energies (not shown) . The domi- 
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Fig. 2. Total excitation energy of the dipoles in the substrate 
Ar384, after deposition of various metal species, as indicated, 
with initial kinetic energy Eq = 0.136 eV/ion, as a function 
of time. 

nant effect is obviously due to charge. In both mass 
cases (Na vs Nag), the dipole energies associated to 
charged projectiles are several orders of magnitude 
larger than for the neutral cases. For neutral sys- 
tems, one can spot a tiny mass effect by comparing 
Na to Nag but which, however, remains negligible 
with respect to charge effects. One also can notice 
a difference between anions and cations, the differ- 
ence lying within one order of magnitude. The neg- 
atively charged electron cloud of the anions experi- 
ence a stronger Pauli repulsion from the Ar cores. As 
a consequence, these projectiles cannot transfer as 
much energy to the Ar substrate as in the positively 
charged cases. This confirms earlier calculations on 
similar systems |36j . 

4. Details of energetics 

From the previous section, we learnt that the ma- 
trix is qualitatively excited the same way either by 
Na + or Na^j". To get more insight into the energet- 
ics, we present in this section the detail of the energy 
sharing for the case of deposition of Na + , which kicks 



out the Ar atom just below the impact point and fi- 
nally stays between the first and the second layers 
of the surface [36] . This case has the advantage that 
no valence electron contributes and that only one 
Na + ion is involved. Thus, there remain only five 
components to the total energy E to t , namely 

J? pmat , 771 , 771 1 7T'Na + , nimat f r,\ 

^tot — -C-pot +-C'dip+A:oupl + -C'ki n +-kkm ■ \ l ) 

The kinetic energies are obvious. The contributions 
to the potential energy are : E™ |f as the potential en- 
ergy of the Ar matrix, consisting out of the Coulomb 
energy and a contribution from the Ar- Ar core repul- 
sion, -Edip as the dipole excitation energy defined in 
Eq. (JT|) , and i^coupi for the coupling energy between 
the Na + and the Ar matrix. As it should be, the 
total energy is well conserved during the deposition 
process (up to a relative error less than 0.01% over 
the whole simulation). Fig. k3] displays various terms 
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Fig. 3. Time evolution of the various contributions to the 
total energy as given in eq. (|2]l in the deposition of Na+ on 
Ar384, with respect to the total initial energy. 

of Eq. (|J as a function of time. At initial time, the 
total energy E to t of course consists of the initial Na + 
kinetic energy but mostly of the matrix potential en- 
ergy -E™t' (99.8 %) with small additional contribu- 
tions from the coupling and dipole energies. These 
two latter contributions do not exactly vanish at ini- 
tial time due to the initial finite distance at which 
the Na + ion is placed with respect to the Ar sur- 
face. Since the initial potential energy surface scales 
with the matrix size (hence with the finiteness of the 
representation of the surface), it makes little sense 
to keep it in the picture. We should hence remove it 
from the total energy (it is anyway a constant). But 
for sake of readibilty of Fig. [3] it appears simpler to 
remove the total initial energy E tot (t = 0) so that 
the components of the energy now sum up to zero at 
any time, as is clear form the figure. All contribu- 
tions have a similar time structure. There is a fast 
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change within the first ps and then the values stabi- 
lize with some final fluctuations. The final share is 
the following : The largest and attractive contribu- 
tion comes from the coupling energy. This is coun- 
terwcighted to comparable parts between matrix po- 
tential energy, matrix kinetic energy, and dipole en- 
ergy (still with the latter taking the lead) . The ki- 
netic energy of the deposited Na + ion is negligible 
as it should be for a well bound particle. 

It is also interesting to show the share of ener- 
gies in relative units, we compare it with the total 
attachment energy of Na + to the surface which is 
-^deposit = —4.7 eV. Relative to the absolute value 
of that energy, we have the contributions : 

77imat 171 TP rp Z^Ar Z?Na 

-^pot ~~ ^tot -C/coupl -C'dip -^kin ^kin 

64% 56% 21.9% 13.7% 0.055% 

The main contributions thus come from the poten- 
tial energy due the rearrangement of the whole ma- 
trix, from the strong coupling between the Na + and 
the matrix (since the metal ion finally locates be- 
tween the two first Ar layers) , and from the internal 
excitation of each Ar atom. 

5. Localization of the dipoles 

The role of charge in the dipole energy suggests to 
explore in more detail the spatial distribution of the 
dipole excitation, with the intuition that it might 
resemble the response to a mere charge. We first 
present in Fig. [4] the dipole energies as a function of 
the axial coordinate p = y/x + y of the Ar cores, at 
impact time and only in the upper layer of the Ar 
substrate, for the case of the deposition of Na + with 
initial kinetic energy of 0.136 eV. The impact point is 
located at p — 0. We observe a high excitation of the 
dipoles which is strongly located around the impact 
point. As we shall see below, one does not observe 
any sizable evolution of this distribution, at least up 
to the times computed here, in terms of total dipole 
energies (see Figs. [T] and 111. In order to analyze this 
aspect in more detail, actual dipole moments (in 

fact, r.m.s. dipole moments d = Jd\ + dy + <P Z ) as 

a function of time are plotted in Figs [5] (case of Na 
monomers) and [6] (case of Na6 clusters). We discuss 
both cases simultaneously because they deliver very 
similar messages and because we have seen that mass 
effect is expectedly small. And yet, it is interesting 
to countercheck by comparing both, since the finite 
extension of Nag might influence the spatial distri- 
bution of Ar dipole excitations. In both figures are 
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Fig. 4. Axial distribution of the dipole energies in the first 
layer of the Ar(001) substrate at impact time, for the depo- 
sition of Na + with initial kinetic energy Eq = 0.136 eV. The 
insert shows a top view of the distribution in the first layer 
(x — y plane). 

plotted four snapshots corresponding to initial time, 
time of first impact and two later instants. There is 
obviously no significant time evolution of the distri- 
butions. One rather observes in Figs. [5] and [6] very 
similar pattern in both late times and both systems. 
One can spot differences in the neutral case; how- 
ever the corresponding values are very small and one 
can probably ignore the point. Initial distributions 
as well are very similar in both figures, at least for 
charged clusters/ atoms. The very regular trend is 
typical of a population of dipoles subject to a dis- 
tant point charge, as is the case for initial states. The 
most interesting panels are probably the ones corre- 
sponding to impact time : One can note a significant 
perturbation of the distributions. This is especially 
true in the case of Na6, most probably because of 
its finite extension which perturbs the matrix on a 
larger range. The effect actually remains for longer 
times, with a somewhat (although not significantly) 
fuzzier distribution in the Nag case. The predomi- 
nant effect is that finally the dipole excitations re- 
main strongly located close to the cluster impact 
point. Everything thus looks as if the charge was 
simply put closer to the surface with the ensuing 
enhanced response of the dipole accompanied by a 
bit of "noise" at the distant points. In that respect, 
the details of the deposition dynamics are to a large 
extent irrelevant. Only charge and its localization 
are really important. This is in accordance with our 
previous study of neutral or charged Na monomer 
deposition/reflection on an Ar surface |36j . Indeed, 
as soon as the Na projectile is reflected, no sizable 
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Fig. 5. Root mean square dipole moments of Ar atoms as function of Ar atomic radius, at four different times as indicated, for 
deposition with initial kinetic energy Eq = 0.136 eV of Na, Na+ and Na — on Ar(001). 



dipole energy beyond "noise" is left in the Ar sub- 
strate. Substantial dipole response is seen only at 
impact time, that is, when the (charged) projectile 
is sufficiently close to the Ar atoms. Thus, in these 
deposition processes, we are facing predominantly 
an electrostatic effect of the polarization of the sur- 
face. This effect is, however, energetically important 
and needs to be taken into account. 



6. An example from embedded clusters 

We have focused up to now on the case of cluster 
deposition. It is also interesting to study what occurs 
in the embedded case in which high charge states 
can be easily attained by laser excitation. We take 
as an example the irradiation of a Nas cluster em- 
bedded inside a finite Ar434 matrix. The laser pulse 
is kept short to concentrate the ionization to a fairly 
well defined initial time and the intensity is tuned 
such that the irradiation leaves the cluster with a 
net 3+ charge. A Coulomb explosion of the clus- 
ter is hindered by the matrix which stabilizes that 
high charge state, but allows a sizable oblate expan- 
sion of the cluster |33l34j . The point is illustrated in 



the upper panel of Fig. [7] which shows the ionic and 
atomic positions along laser polarization axis as a 
function of time. The Ar matrix is also strongly per- 
turbed with significant atomic rearrangement but 
the whole system finally remains stable. The lower 
panel of Fig.[7]displays energies : The kinetic energy 
of Na ions and Ar cores, and the energy stored in 
the Ar dipoles. The energy balance is to some ex- 
tent incomplete as it misses the potential energy of 
cluster and matrix. But, as in case of deposition, the 
kinetic and dipole energies provide sufficient infor- 
mation about energy flow. We see in the lower part 
of Fig. [7] a behavior similar to the deposition case, 
namely that the initial Na motion is quickly damped 
and its energy accordingly transferred to the matrix. 
Again, one observes a significant amount of energy 
stored in the Ar dipoles, about three times as much 
as in the Ar core kinetic energy. 

To complement the analysis of this irradiation 
case, we consider again the spatial distributions of 
dipoles at various instants in Fig. [8] The geome- 
try is now a bit more complicated and we should 
care for both directions, along the laser polariza- 
tion (z direction) and perpendicular (axial direc- 
tion) to it. We thus slightly change the representa- 
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Fig. 6. Root mean square dipole moments of Ar atoms as function of Ar atomic radius, at four different times as indicated, for 
deposition with initial kinetic energy E = 0.136 eV/ion of Na 6 , Na+ and Na~ on Ar(001). 



tion of the distribution in Fig. [8] considering sepa- 
rately both directions. The situation is also differ- 
ent from the charged deposition case to the extent 
that the cluster is initially neutral. This explains the 
rather "democratic" distribution of dipoles all over 
the matrix at initial time (upper panels). This con- 
firms earlier findings on the role of Ar polarizabili- 
ties even in low energy phenomena such as optical 
response }29|32| . The initial dipole amplitudes are 
rather small as compared to the values which they 
acquire when the cluster becomes charged (mind the 
scales of upper and lower panels). The lower panels 
of Fig. [8] show large values of dipoles, comparable 
to deposition of charge clusters (see Figs. [5] and [6]). 
The distribution, though, looks at first glance quite 
different from the deposition cases with a maximum 
at a finite distance from center rather (in radial co- 
ordinates, right lower panel) than a monotonously 
decreasing shape. The effect is in fact due to the fi- 
nite (large) extension of the Na cluster after irradi- 
ation. Indeed, a sizable fraction of Ar sites are em- 
bedded in the Na cluster electron cloud and thus see 
a screened charge, whence the reduced dipole po- 
larization. In order to exemplify the point, we have 
also plotted in Fig. [8] the actual position of the most 



external Na ions and this makes the pattern clear, 
with an effect only along radial coordinate due to the 
oblate shape of the irradiated cluster |33|34j . Apart 
from that detail, we find again a dominantly static 
charge effect explaining the dipole polarizations of 
Ar atoms, much similar to the deposition case with 
localization prevailing again and little own dynam- 
ics. 



7. Conclusion 

We have discussed in this paper the atomic re- 
sponse of an Ar surface or matrix perturbed by a 
metal cluster cither through deposition or irradia- 
tion by a laser. We have focused the analysis on 
the internal response of the Ar atoms by studying 
their single internal degree of freedom in our model, 
namely their dipole polarizability. We have seen that 
the dominant effect is due to cluster charge. Mass 
effects are mostly negligible with respect to charge 
effects. We have also analyzed how the dipole exci- 
tations are spatially distributed. We have observed 
that the excitation remains strongly localized and 
essentially does not evolve in time. It resembles in 
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Fig. 8. Root mean square Ar dipoles for the hindered Coulomb explosion of Nag embedded in Ar434, exposed to a laser of 
intensity 2 X 10 12 W/cm 2 , frequency u) = 1.9 eV, and FWHM=33 fs. Left panel : Distribution as a function of the Ar z 
coordinates; right panels : That as a function of the Ar axial coordinates p = \/ x 2 + y 2 . Top panels : At initial time; bottom 
panels : For three subsequent times as indicated. The maximum of excitation energy observed in the bottom right panel at 11 
ag is due to the oblate deformation of the created Nag + . The vertical lines in the bottom panels indicate the corresponding 
coordinates of the outer Nag ions. 



many respect the effect of a finite charge deposited at 
some place in the system. This holds true, up to de- 
tails, for all situations involving charges, in deposi- 
tion as well as in irradiation dynamics. Independent 
dynamical evolution of the Ar dipoles creates some 
background noise which, however, remains quanti- 
tatively unimportant at the present level of analy- 
sis. The quantitatively very important effect remains 
the nearly static dipole deformation. The Ar dipoles 
store a sizable fraction of the available energy. This 
is in perfect qualitative agreement with recent ex- 
periments on similar systems. 
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